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ABSTRACT 

X-ray sources in spiral galaxies can be approximately classified into bulge 
and disk populations. The bulge (or hard) sources have X-ray colors which are 
consistent with Low Mass X-ray Binaries (LMXB) but the the disk sources have 
softer colors suggesting a different type of source. In this paper, we further 
study the properties of hard and soft sources by constructing color segregated X- 
ray Luminosity Functions (XLF) for these two populations. Since the number of 
sources in any given galaxy is small, we coadded sources from a sample of nearby, 
face-on spiral galaxies observed by Chandra ge Project in Cycle 2. We 

use simulations to carefully correct the XLF for completeness. The composite 
hard source XLF is not consistent with a single power-law fit. At luminosities 
Lx > 3x10^^ erg s~^ it is well fit by a power law with a slope that isconsistent 



with that found for sources in elliptical galaxies by lKim &: Fabbiand (l2004l ). This 
is supports the suggestion that the hard sources are dominated by LMXBs. In 
contrast, the high luminosity XLF of soft sources has a slope similar to the 
"universal" High Mass X-ray Binary XLF. Some of these sources are stellar mass 
black- hole binaries accreting at high rates in a thermal/steep power law state. 
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The softest sources have inferred disk temperatures that are considerably lower 
than found in galactic black holes binaries. These sources are not well understood, 
but some ma y be super-soft u ltra-l uminous X-ray sources in a quiescent state as 
suggested by lSoria fc Ghosh 



1. Introduction 



It has been known for three decades tha t nearby spiral galaxies contain a multitude 
of X-ray sources (Ivan Speybroeck et al.lll979l ). Observations of the Milky Way and Local 
Group galaxies have shown that the brightest of these sources are Low Mass X-ray Binaries 
(LMXB). In these systems, material is transfered at a high rate via Roche Lobe overflow 
from a low mass star onto a compact companion (a white dwarf, neutron star or black hole). 
Their X-ray s pectra are broadly characterized by a power law with photon index 7 ~2.5 
(iTanakalllQQTI ). High Mass X-ray Binaries (HMXB) are powered by Bondi-Hoyle accretion 
from the stellar winds of a young, high mass star onto a compa ct object. They typicall y 
have much harder spectra than LMXB (photon index 7=0-1.5, iHaberl &: Pietschl (120041 )) 
The accretion rates are lower than for Ro che Lobe overflow sources and the corresponding 
X-ray luminosities also lower (jWhitelll989l ). In addition to the basic classification depending 
on the mass of the secondary, many other groups of sources have been identified depending 
on their X-ray spectra and variability properties. For example, so-called Super Soft Sources 
(SSS), believed t o be binaries with white dwarf prima ries accreting matter at a highly super- 
Eddington rate (IGreiner. Hasinger. fc Kahabkalll991l) have essentially all of the ir flux below 



IkeV, and are highly variable ( Di Stefano fc KongI 2003 : Di Stefano et al. 2004). Supernova 



remnants are also kno wn to be s t rong X -ray sources (iReynoldd l2008l ) . For a review of X-ray 
sources in galaxies see iFabbianol (120061 ) . 



In the closest galaxies, clues as to the n ature of X-ray sources can be derived from th e 
optical environments and/or counterparts (IKilgard et al.l l2005l : iTrudolyubov et al.l l2005l ). 
However, it is more difficult to identify X-ray sources in galaxies which are too distant for 
optical counte rparts to be dete c ted a nd where there are not enough counts to get an X-ray 
spectrum. In iPrestwich et al.l (120031 ) we suggested that X-ray colors are a good starting 
point for source classification. We found that the X-ray colors of the sources described 
above (low and high mass binaries, SSS, and SNR) lie in different parts of the X-ray color- 
color diagram, albeit with some overlap of source types. In particular, sources in elliptical 
galaxies and spiral bulges fall in a well defined region of the X-ray color-color diagram. We 
suggest that this population is dominated by Low Mass X-ray Binaries. We find very few 
sources with hard HMXB colors, which is to be expected given their low X-ray luminosities. 
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We also find a population of faint X-ray soft sources in disk galaxies not seen in bulges. 
These soft/thermal sources are clearly associated with star for mation, and are l i kely a mix 
of supernova remnants and accretion powered sources (see also Di Stefano et al. (j2004 )). 



In this paper, we further investigate the properties of hard (or LMXB) and soft sources 
by constructing color segregated X-ray Luminosity Functions (XLFs). Our motivation is to 
confirm our interpretation that the hard sources are dominated by LMXBs and to better 
understand the mix of objects that fall into the soft class. In Section [2] we describe our 
galaxy sample and completeness simulations and in Section [3] we describe fits to hard and 
soft sources. In Section H] we compare our results to previous work and discuss the origin of 
soft sources. Our conclusions are summarized in Section [61 



2. Galaxy Sample and Simulations 



Most spiral galaxies have 10-50 X-ray sources, so that it is impossible to construct 
statistically significant color segregated luminosity functions for a single galaxy. We therefore 
constructed color-segregated XLFs by coadding hard and soft sources from 11 spiral galaxies. 
The galaxies chosen for this work were observed as a Chandra Large Project in Cycle 2 (see 
Kilgard et al.l (l2005l ) for details of data reduction and a complete source list). They are 
all nearby (< 10 Mpc), face on spirals with low foreground absorption. The X-ray colors 
of sources in this sample are therefore relatively insensitive to patchy local absorption and 
the inclinatio n of the galaxy. We use color definitions hard color (HC) and soft color (SC) 
introduced by lPrestwich et al.l (120031 ): 



HC 



SC 



H-M 
T 

M-S 
T 



Here if, M, and 5* are counts in the soft (0.3-1 keV), medium (1-2 keV) and hard (2-8 
keV) bands respectively and T = S + M + H . We define soft and hard source colors as 
follows: 

soft color: SC <QA 

hard color: -0.4 < SC < 0.2, -0.3 < HC < 0.4 
These regions are shown on a color-color diagram in Figure [H 

The composite XLFs for soft and hard sources sources are shown in Figures [2l and 
|3l Inspection of these figures suggest that the soft sources are well represented by a single 
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power law. In contrast, the XLF of the hard sources starts to flatten at ~ 3x10^^ erg s~^. 
The flattening of the XLF may be due to a real change in the properties of the sources or 
an artifact due to incompleteness (or a combination of both). There are several factors that 
affect the detection threshold for X-ray sources, including the presence of diffuse emission 
and the distance of the source from the aimpoint. In addition, ACIS-S was significantly more 
sensitive to soft emission (< 1 keV) in Cycle 2 than to harder sources; this will result in a 
different completeness limit for hard and soft sources. We construct simulated XLF to help 
distinguish between a real break in the XLF and incompletemess effects. 



We use the "backward" method described by iKim fc Fabbiand (l2004j ) to ascertain the 
completeness limit of each galaxy for a given source type (hard or soft). This method involves 
simulating a source with the appropriate spectrum, and adding the simulated source to the 
real galaxy data. The source detection algorithm used for the original galaxy source list is 
then run again to determine whether the "fake" source is detected. Approximately 20,000 
soft sources and 20,000 hard sources are added to each galaxy, sprinkled randomly over the 
area within D25 . We use a power law of photon index 1.2 to simulate hard sources and a 
power law of photon index 3.0 to simulate soft sources. The detect algorithm was run on 
the full band image, as was also the case for the galaxy analysis. The ratio of the number of 
input to the number of detected sources at a given luminosity is the key factor: when this 
falls below 100% the sample is not complete. In order to construct a composite XLF, we 
added sources to luminosity bins where the completness fell below 100%. For each galaxy, 
sources were added to the XLF according to the following formula: 



NaM{L) = ^-^N^L) (1) 

Here Nadd{L) is the number of sources added to the bin of luminosity L and F is the 
fractional completeness, defined as the ratio of the number of input to the number of sources 
detected in simulations. Ngai{L) is the number of real sources in the galaxy with luminosity 
L. 

We find that the soft source population is 100% complete down to 3.9x10^^ erg s~^ 
for most galaxies (MlOl, IC 5332, M94, M51, M83, NGC 2681, M74 and NGC 1291). No 
sources were added to the XLF for these galaxies. We omit galaxies NGC 4314, NGC 3184 
and NGC 278 from the soft source composite XLF because they are less than 40% complete 
at 3.9x10^'' erg s"^. Detection of hard sources was complete down to 7.5x10^"^ erg s~^ for 
galaxies IC 5332, MlOl, M51, M83, NGC 2681, NGC 278 and NGC 628. Sources were added 
to the lowest luminosity bins of the remaining galaxies as summarized in Tabled! In addition 
to correcting for completeness in the lowest luminosity bins, we also removed bright nuclear 
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sources from the complosite XLFs. In many cases these are weak LINERS and not X-ray 
binaries. 



3. Fits to Composite XLF 



The method used for fitting is described by iKilgard et al.l (120051 ) and full details can be 
found in that paper. Here we give a brief summary of the procedure. 

We fit a single p ower law to the unbinn ed, differential, composite XLFs using a maximum 
likelihood statistic (ICrawford et al.lll970l ). The XLFs were corrected for completeness as 
described in the previous section. The single power law has the functional form: 



dN 



dL 



X 



r-7 



(2) 



The goodness-of-fit estimate is performed by simulating a luminosity distribution with 
the best fit slope. One million iterations were performed. If the data are well fit by a single 
power law the goodness-of-fit statistic (GOF) approaches 1.0. The contrib ution of cosmic 



backg round sources to the XLF was evaluated using the LogN-LogS curves of iGiacconi et al. 
( 120011 ). The soft band was used for the soft sources and the hard band for the hard sources. 



We fit down to the completeness limit (3.9x10^^ erg s~^) for soft sources. It is impossible 
to fit a single power law from the completeness limit to the highest luminosity source for 
the hard sources. We therefore fit a single power law above 3x10^^ erg s~^. Figures [2] and 
[3] show the composite XLF and best fit power law (black curves) for soft and hard sources. 
Also shown is the slope representing the background source distribution (red curve) and the 
background subtracted data and best fit (blue curve). Best fit parameters to the composite 
XLF are shown in Table [21 



4. The X-ray Source Population in Spiral Galaxies 

There are two basic conclusions to be drawn from the results presented in the previous 
section. The first is that although the completeness limit for a given observation depends on 
many factors, such as diffuse emission, it is lower for soft sources than for hard sources. The 
second is that there appears to be a break in the hard source XLF at about ~3xl0^^ erg 
s~^. The flattening of the hard source XLF shown in Figure [3] at low luminosities (<7xl0^^ 
erg s~^) is largely due to incompleteness. 
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4.1. Hard Sources 



In iPrestwich et al.l (120031 ) we suggest that the hard sources are dominated by Low Mass 
X-ray Binaries. The XLF of sources in elhptica l galax ies, which are a lmost e xclusively LMXB 
has be en studied in depth by lKim fc Fabbiand (120041 ) and lGilfanovl (120041 ) . iKim fc Fabbiano 
( 120041 ) find that the XLF of a sample of 14 elliptical galaxies (corrected for incompleteness) 
can be fit with a single power law, 7=1.9 over the range 2x10^"^ — 2 x 10^^ ergs s~^. Although 
not formally required, a broken power law gives a better fit (71=1. 8±0. 2, 72=2. 8±0. 6) with a 
break luminosity Lb=5x 10^^ ergs s~^. Our fits to the hard sources show very similar features: 
an XLF which steepens at a few x 10'^'^ ergs s^^. Our fit parameters are consistent with those 
of iKim fc Fabbiand (120041 ). albeit with large u ncertainties due to a sma ller sample. It is likely 
that we are sampling the same population as iKim &: Fabbiand (120041 ). 



4.2. Soft Sources 



In IPrestwich et al.l ( l2003l ) we suggested that many of the soft sources are thermal 
supernova remnants ( SNR). We stress here that most of the soft sources considered in 
Prestwich et al.l (120031 ) have luminosities <1.0xlO^^ erg s~^. They are below the complete- 
ness limit for the composite XLF derived here. The soft sources at the high luminosity end 
of the XLF are almost certainly dominated by accretion powered sources. The association of 
soft sources with star formation strongly suggests that the donor stars are young supergiant, 
O or B stars, even though the soft colors of these binaries is very different from the hard 
spectra of wind-accretion sources described in Section [H We therefore compare the XLF of 
soft sources to that of HMXB. The nature of the soft sources is discussed further in Section |5] 

The XLF of high mass b i naries in star forming galaxies has been studied in detail by 
Grimm. Gilfanov. fc SunyaevI (120031 ). They find that the HMXB XLF is well described by a 
single pow er law (7 = 1.61±0.12) over the lurn inosity range 10^^ — 10^° ergs s~^. The slope 
derived by iGrimm. Gilfanov. &: SunyaevI (120031 ) is very close to our best fit single power law 
for soft sources (7 = 1.73 ± 0.15), suggesting that we are sampling the same population. 



5. The Nature of the Soft Sources 

As discussed above, the luminous soft sources discussed here are most likely HMXB. 
However, they are clearly more luminous and have much softer spectra than wind-fed HMXB 
(binary pulsars) found in the Milky Way and local group galaxies. The higher luminosities 
suggest that the soft sources are accreting material at a higher rate than is found in wind-fed 
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HMXB. This requirement, plus the softer spectra, suggests that the accretion mechanism is 
fundamentally different than galactic HMXB. In this section we discuss the nature of these 
sources in more detail. 

The association of the soft sources with star formation is clearly demonstrat ed in the case 



of NG C 4736 (M94). This galaxy has a circumnuclear ring of star formation (jWaller et al. 



( l200ll )). Figure m shows the Ha image of this galaxy, with the positions of soft (green) and 
hard (blue) X-ray sources. It is clear that the soft sources are found primarily in the star 
forming ring. Given the association with star formation, these sources are most plausibly 
neutron star or black hole X-ray binaries with young stellar compar iions. The X-ray spectral 



states of black hole X-ray binaries have been described in detail by lRemillard &: McClintock 



( 120061 ). In the low accretion state the X-ray emission is dominated by a hot corona. As 
the accretion rate increases the X-ray flux comes primarily from a disk and is thermal in 
origin. The emission is characterized by a disk-blackbody spectrum (kT~l keV) with some 
contribution at higher energies from a power law component. In the highest accretion state 
(Steep Power Law, or SPL) the X-ray emission is characterized by a soft power law (P >2.5). 
Figure [5] shows the X-ray color-color diagram from ?? for all detected sources, with thermal 
and SPL models plotted for reference. The red curve shows the colors of disk-blackbody 
models for inner disk temperatures of 0.1-1.1 keV. Adding a power law component to a 
disk-blackbody has the effect of moving the colors up and to the left of the diagram, as 
demonstrated by the green points. The purple points show a simple power law, representing 
the SPL state. 

Figure [5] shows that some of the soft sources have colors that are very similar to black 
holes accreting in the thermal or SPL state. However, some of the softest have disk temper- 
atures that are much cooler (kT~0. 1-0.3 keV) than is typically found in galactic black hole 
binaries (kT~l keV). Low disk temperatures have also been observed in Ultra-Luminous X- 
ray Sources (ULX). Since M ~ T^~^, such low disk temperatur es have been used to sug gest 



that ULX contain Intermediate Mass Black Holes (IMBH, e.g. iMakishima et al.l ((20001)). It 
seems very unlikely that we observing a population of IMBH. If the mass of the black holes 
are ~ 100 Mq , the low luminosities (10^^ — 10^* ergs s~^) would imply sub-Eddington ac- 
cretion rates. In which case, we would not expect to see thermal emission, which is generally 
observed when the accretion is dominated by a disk and the rates are close to Eddington. 
The combination soft spectrum and low luminosity displayed by a subset of our soft sources 
is hard to explain if the sources are IMBH. If they are stellar mass black holes, their disk 
temperatures are unusually low. 



Soria &: GhoshI (120091 ) recently discussed a state transition in a ULX in NGC 4631. In 



the high state NGC 4631 is a super-soft source. In a low state observed by Chandra it has 
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a luminosity ~ 10^^ — 10^^ ergs s~^ and colors consistent with a thermal temperature of 
0.1-0.3 keV. They interpret the ULX/super-soft phase of NGC 4631 as a super-Eddingon 
outburst powered by nuclear burning on the surface of a white dwarf. In it's low state, 
the NGC 4631 ULX is very similar to the soft sources described here. By extension, some 
soft sources in spiral galaxies may be quiescent super-soft ULX. A detailed study of the 
variability properties of these sources is required to better understand their nature. 



6. Summary and Conclusions 

In this paper, we construct color segregated X-ray Luminosity Functions (XLF) of hard 
and soft sources by combining data from a sample of face-on spiral galaxies. Since the number 
of sources in any given galaxy is small, we coadded sources from a sample of nearby, face-on 
spiral galaxies. To determine the compleness level of each galaxy as a function of source 
type (hard or soft) we add "fake" sources to real galaxy data and determine the ratio of fake 
sources added to the number recovered by source detection routines. If all fake sources were 
detected the galaxy is 100% at that flux level. We added sources to the lowest luminosity 
bins to produce composite color segregated XLF that are complete down to 3.9x10'^^ erg s~^ 
for soft sources and 7.5x10^^ erg s~^ for hard sources. 

The composite hard source XLF is best fit by a broken power law an d has a high lumi- 



nosity slope consistent with that found for sources in elliptical galaxies by lKim fc Fabbiano 



(120041 ). We conclude the hard population is dominated by LMXB. The soft XLF above the 



compl eteness limit is very similar to the "universal" HMXB XLF described by lGrimm. Gilfanov. fc Sunyaev 



(120031 ). This strongly suggests that although thermal supernova remnants contribute to the 
soft population at lower luminosities, the high luminosity end of the XLF is dominated by 
accretion-powered HMXB. The association of soft sources with star formation is confirmed 
by looking in detail at NGC 4736, where the soft sources are coincident with a ring of 
Ha emission. Some of these sources are stellar-mass black hole binaries accreting in the 
thermal/PL states. However, the softest sources have inferred disk temperatures that are 
considerably lower than those found in galactic black holes binaries. These sources are not 
well understood, but some may be super-s oft ultra-luminous X-ray sources in a quiescent 



state as suggested by lSoria fc GhoshI (120091 ). 
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Fig. 1. — X-ray color-color diagram for sources in spiral galaxies from ??. The location of 
hard and soft sources is shown 
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XLF for Soft Sources 



100 



A 10 



1 



~i 1 n 1 — I — I — TT 



1 1 1 1 — I — I — n~ 



~i 1 1 1 — I — I — r~r 




10 



37 



1 0=^' 1 0^' 
Luminosity [erg/s] 



10 



40 



Fig. 2. — XLF for soft sources. The black curves show the data and best fit power law, the 
red line shows the distribution of background sources and the blue curves the background 
corrected XLF. The vertical dotted line shows the completeness limit. Points below the 
completeness limit were not included in the fit. 
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XLF for Hard Sources 




Luminosity [erg/s] 

Fig. 3. — XLF for hard sources. The black curves show the data and best fit power law, the 
red line shows the distribution of background sources and the blue curves the background 
corrected XLF. The red vertical dotted line shows the completeness limit. The black vertical 
line shows the lowest luminosity included in the fit. 
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Fig. 4. — Ha image of the center of NGC 4736 showing the star forming ring. Positions of 
soft X-ray sources are plotted in green and hard sources plotted in blue. The soft sources 
are predominantly associated with the ring. 
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Fig. 5. — X-ray color-color diagram for sources in spiral galaxies. Curves for thermal and 
SPL models for black binaries are plotted for reference. The red curve shows the colors 
of disk-blackbody models for inner disk temperatures of 0.1-1.1 keV. Adding a power law 
component to a disk-blackbody has the effect of moving the colors up and to the left of the 
diagram, as demonstrated by the green points. The purple points show a simple power law 
characteristic of sources in the SPL state.. 
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Table 1. Completness Limits for Hard Sources 



Luminosity 


NGC 4314 


NGC 1291 


NGC 3184 


M94 


xlO^^ erg s~^ 










0.75 


20% (5) 


50% (4) 


50% (1) 


70% (1) 


0.95 


50% (2) 


68% (1) 


70% (1) 


90% (1) 


1.15 


60% (1) 


77% (1) 


80% (0) 


100% (0) 


1.35 


80% (0) 


90% (0) 


90% (1) 


100% (0) 


1.5 


90% (1) 


100% (0) 


100% (0) 


100% (0) 



Note. — Table gives completeness as a percentage of input sources 
detected. The number in parentheses is the number of sources added 
to the XLF to account for completness. Galaxies not hsted are com- 
plete down to 7.5x10^''' erg s"-^. 



Table 2. Power law fits to composite, corrected XLF 



Sample 


7 


'Ycorr 


lO^'' erg 


lO^'' erg 


GOF 






Ntot 


Soft 


1.82±0.14 


1.73±0.15 


3.9 


3.9 


0.7 


32 


32 


170 


Hard 


2.53±0.25 


2.55±0.28 


30.0 


7.5 


1.05 


37 


104 


236 



Note. — Power law fits to composite, corrected XLF. 7 is the fitted slope to the XLF with no background 
correction, jcorr correcting for background sources. Li^^ is the low luminosity used in the fit, Lcomp the 
completeness limit. GOF is the goodness-of-fit parameter and Nfn the number of sources included in the 
fit. Nc omp is the number of sources above the completeness hmit and Ntot is the total number of detected 
sources. 



